General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



FINAL TECHNICAL REPORT 


entitled 


Atomic Oxygen Between 80 and 120 KM 
Evidence for a Rapid Spatial Variation 
In Vertical Transport Near the Ionosphere 


(MASA>CB-155€09) AICBIC ClXGfM fiFTHlEN 80 
AND 120 ka: EVIDENCE FOE A RAPID SPAIIAL 

VABIAIION IN VEEIlCAl TRANSPORT NEAR THE 
IONOSPHERE Eioal Tccfaoical Report* Sep. 
1976 - Auq. 1977 fHicfaiqan Ooiv.l 46 p HC 


N77-32660 

fit 

One las 

G3/46 47690 


prepared under 


NASA Grant NSG-7279 

Title: Analysis of OGO-6 Photometer Data 
Period: September 1976 - August 1977 




by 


Thomas M. Donahue 
Barry Wass*»r 

Department of Atmospheric and Oceanic Science 
Space Physics Research Laboratory 
The University of Michigan 
Ann Arbor, Michigan 48109 





ATOMIC OXYGEN BETWEEN 80 AND 120 KM: EVIDENCE FOR A RAPIL 

SPATIAL VARIATION IN VERTICAL TRANSPORT NEAR THE MESOPAUSE 


Barry Wasser and T. M. Donahue 

Department of Atmospheric nnd Oceanic Science 
Space Physics Research Laboratory 
The University of Michigan 
Ann Arbor, Michigan 48109 


ABSTRACT 


Analysis of OGO-6 01 Green Line photometer results has been carried out 
for 8 cases when the alignment of the spacecraft was such that local emission 
rates could be determined below the altitude of maximum emission and down to 
about 80 km. The results show a variation on a scale of 6" - 8* in latitude 
between regions where the emission rate increases rapidly between 90 and 
95 km and regions where it Increases slowly from 80 km to 95 km. Latitude- 
altitude maps of iso-emlsslvlty contours and Iso-density contours for [0] are 
presented. The latter are computed under three assumptions concerning excita- 
tion mechanisms. Comparisons of the spatial variations of oxygen density with 
the results of a time dependent theory suggest the regions of strong downward 
transport alternate on a scale of about 1000 km with regions of weak transport 
near 90 km. In the first case conversion of 0 to 0^ at night appears to be 
overwhelmed by downward transport of 0. 


INTRODUCTION 


Slant column emission rates from the atomic oxygen green alrglow at 557.7 
nm were measured by the horizon scanning photometer flown on OGO-6. Reports 
of results obtained from the analysis of these observations have been published 
In a series of papers [ Thomas and Donahue , 1972; Donahue et al ., 1973, 1974; 
Donahue , 1975; Donahue and Carignan , 1975]. A detailed description of the 
instrument and its operation can be found in Donahue et al . [1973] . The 
photometer was equipped with a stepping mirror so that it could scan the hori- 
zon in 128 steps of 7.5* from 10^ to 26° below the local horizontal and had an 
angular field of view of 7.5' in the vertical and 4.2® horizontally. The 
satellite orbit had a perigee of 397 km and an apogee of 1098 km and an in- 
clination of 82°. From perigee the field of view of the photometer encompassed 
a vertical surface less than 5 km in extent 100 km above the earth's surface. 

The green alrglow features normally detected by the photometer were an 
region glow at tropical latitudes associated with the electron recombina- 
tion of 02 "*" [ Thomas and Donahue , 1972] and the well known alrglow layer in 
the lower thermosphere generally assumed to be associated with the recombi- 
nation of atomic oxygen. The excitation mechanisms proposed for this alrglow 
are usually associated with the names of Chapman [1931], Barth [1964], Slanger 
and Black [1977] 


or 


0 + 0 + 0->02 + 0(‘S) 
0 + 0 + M->02*+M 

0 + O 2 * O 2 + 0(*S) 


( 1 ) 

( 2 ) 

<3) 


followed after 0.75 seconds by 


0(^S) 0(’P) + hv 


(A) 


Alternatively, 0(^S) can be quenched by two processes 


0(*S) + 0(^P) 0(3p) + 0(^P) (5) 

0(^P) + 0(*D) 


or 


0(^S) +02-^ 0(*P) + O 2 (6) 

In the early analysis of the photometer data [ Donahue et al . , 1973] slant 
emission rates were converted to local volume emission rates by a simple In-* 
version which was based on the assumption that the emission rates at a given 
altitude were independent of latitude. This assumption waA known to be In- 
correct, for the results of the analyses themselves demonstrated large varia- 
tions in the emission rates with latitude. It was characteristic of the ob- 
servations that fluctuations of the order of 50 to lOOZ occurred on a latitude 
scale of tens of degrees superimposed on systematic patterns hemispherical In 
scale. The "uniform layer" assumption nevertheless is a reasonable approxima- 
tion to use in deriving local emission rates from the Integrated slant rates 
for the region above the altitude of maximum emission because the overwhelming 
contribution to the slant emission there comes from a short horizontal region 
near the lowest altitude attained along the line of sight. For fields of view 
that pass below the maximum, however, there are large contributions to the 
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slant intensity from regions closer to and further from the satellite than the 
lowest altitude along the line of sight. This situation becomes more and more 
serious as the altitude probed decreases since the volume emission rate decreases 
with altitude. There is no way to correct for this effect as long as the field 
of view lies in a plane that makes a large angle with the orbital plane. 

However » there were occasions when the satellite yawed in such a way that the 
field of view lay near the orbital plane long enough for a large number of ver- 
tical scan of the photometer's stepping mirror to be obtained. Under these 
circumstances, it is possible to use a technique developed by Thomas and 
Donahue [1972] to obtain reliable estimates of the emission rates at altitudes 
below the location of the maximum emitting region. The reason is that each 
vertical scan occurs rapidly enough that a given point in space is sampled 
many times along varying lines of sight as the satellite moves in its orbit 
(Figure 1) . Figures 1 and 2 in the paper by Thomas and Donahue [1972] are 
useful in producing a visualization of this effect. 

The latitude spacing between scans such as 1 and 2 or 3 and 4 (Figure 1) 
is about 3.3* at low latitudes and between scans 2 and 3 is only 1.1*. This 
pattern is the consequence of the fact that pairs of scans were performed 
with a green and a yellow (589 nm) filter alternately in the photometer. 


DATA ANALYSIS 


To calculate the local volume emission rates from the data obtained » the 
region studied by a sequence of vertical photometer scans is divided into a 
matrix vhose dimensions are 1 km in altitude by 0.5* in latitude. In the 
case of each vertical scan, the local volume emission rates are calculated 
for a sequence of matrix elements at the latitude where the line of sight 
along which the slant emission rate is a maximum passes closest to the earth's 
surface. These calculations are made using the uniform layer assumption. Next, 
the elements in each row (latitude) and column (altitude) containing an clement 
with a "measured" emission rate are assigned rates determined by linear inter- 
polation between pairs of "measured" emitting elements. All other elements 
are assigned an emission rate of 1 cm”’s~*. This is the first order approxima- 
tion to the emission rate matrix. 

Next, the calculated slant emission rates from the first order matrix 
are determined from the sum over all elements along the lines of sight corres- 
ponding to given mirror positions. If the calculated values differ from the 
observed ones, the emitting elements Involved in each of such sums are corrected 
by multiplying each of them by the ratio of the observed to the calculated 
slant emission rate for the mirror position involved. (Some of the unit 
matrix elements will be involved in this procedure. Their values will then 
be the ratio of observed to calculated slant emission rates.) This procedure 
is followed for each scan of- the photometer mirror. Where an emission rate 
matrix element already traversed by a previous scan is encountered in a sub- 
sequent scan It is omitted from the summation and correction process in the 
latter scan. A second approximation to the emission matrix is thus obtained. 
This matrix is then smoothed by taking a running average in the vertical (in- 
volving three adjacent elements) and in the horizontal (using five adjacent 
elements). During this smoothing process a large number of the elements in 
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the matrix will be assigned values differing from unity. Next, the slant emis- 
sion rates from the smoothed matrix are calculated, compared to the observed 
emission rates, corrected and smoothed. This procedure is normally repeated 
25 times and the result used to plot iso-emlsslon rate maps such as shown in 
Figures 2a and 9a. In the case of pass 2269 on day 313 (Figures 8a and 10) the 
adequacy of the Iteration procedure was calculated by iterating 300 times in- 
stead of 25 times. While some structure in the coarser map of Figure 8a dis- 
appeared from Figure 10 as a result of 12 times as many iterations, the improve- 
ment is not great enough to warrant the 12 fold increase in computer time in- 
volved. 


' t 


RESULTS 


The emission rate maps in Figures 2a-*9a are for the following orbit segments; 


Day 304, 

pass 

2135 

from 

-26.6* 

to 

- 7.1" 

latitude 

304, 

pass 

2142 

from 

-20.1* 

to 

- 4.8" 

latitude 

311, 

pass 

2238 

from 

19.2" 

to 

33.6" 


311. 

pass 

2239 

from 

15.5" 

to 

M.3" 


311. 

pass 

2243 

from 

14.2" 

to 

20.5" 


311, 

pass 

2244 

from 

8.8" 

to 

22.4* 


311, 

pass 

2269 

from 

18.6" 

to 

27.8" 


329, 

pass 

2505 

from 

48.8" 

to 

54.1* 



From the emission rates the atomic oxygen density profiles may be calcu- 
lated provided all relevant rate constants are known and the excitation of 0(^S) 
occurs as postulated. Unfortunately, there is still some uncertainty with re- 
gard to the values for the rate constants, k^, kjj, fc^, k^, and kg. Donahue et 

al. [1973] in analyzing their data made various assumptions with regard to the 
— — ^ 

values of rate constants ka, k^» and kg. One of these was that these rate con- 
stants have the following values. 

k = 4.8 X lO'^^cm^s-^ 
a 

kj - 1.2 X 10“ exp(-226/T) cm’s“^ 

k - 4.3 X 10"“ exp(-853/T) cm^s"* 
e 

where k, and k are taken from Slanger and Black [1973] and k_ from Felder 

d e — a — 

and Young [1972]. Recently, as a result of analyzing data obtained from an 
experiment in which photometer and mass spectrometer observations obtained by 
Atmosphere Explorer (AE-C) were coordinated with a vertical profile for the 
green line emission rate obtained from a sounding rocket [ Wasser et al ., 

1976] deduced a rate constant for k^ given by 
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kg - 2,4 X 10“’^ exp(10QQ/T) cm*8“^ 


Recently, Slanger and Black [1977] have argued from laboratory measurements and 
analyses of satellite and rocket data that the Barth mechanism can account for 
the alrglow. According to their mechanism 
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The emission rates have been converted to oxygen density profiles for each 
of these three sets of assumptions. The cases will be labelled and SB. 

The lso~density contours corresponding to the eight emission rate maps are 
plotted in Figures 2b, c, d-9b, c, d. Vertical cuts through the emission rate 
and density maps at selected latitude are shown in Figures 11-19. The 
generally smaller value of the Chapman rate constant for Case C 2 than for C^ 
and its decreasing value with T causes the C^ densities to be larger than C^^ 



by a factor of about 2 near 95 kin and the ratio of to densities to Increase 
with altitude as T Increases. At 120 km, the ratio Is typically 5. Thus, 
gives larger values for 10] than and a slower decrease with altitude. The 
SB densities are typically smaller than or of the order of densities near 
90 km but rise much more rapidly than densities so that the ratio Increases 
from less than 1 to a value of 2 to 3 at about 102 km. Thereafter the ratio 
of SB to densities decreases very slowly or remains relatively constant to 
120 km at a value in the neighborhood of z. The low densities for a given 
prediction by the SB mechanism near 80 km and the increase with altitude rela- 
tive to is a consequence of the increasing rate for Reaction (2) as [M] 
increases. In fact it is easy to show that 

[OK. k. 1/3 K 

Toif "<rk> ^ V 1 

\ 

Not all of the emission rates obtained by this method are physically 
meaningful. One important property of the photometer detector is that it was 
designed to produce an output volta'^e that stepped in intervals equivalent to 
lOR. Thus a reading of zero meant that the emission in the field of view was 
less than 5R. The first step in output voltage above zero indicated an emis- 
sion rate between 5 and 15R, the second step 15 and 25R and so on. Furthermore, 
it Is necessary to extrapolate a background counting rate from data obtained 
when the field of view was above 120 km and subtract this background from the 
slant rate measured at low altitudes. The unct’rtainty Introduced by background 
subtraction between 90 and 80 km can hardly be less than 20R except in the very 
best of circumstances. It can never be less than 5R. At high altitude the 
emission rate decreased with a scale height typlcaly of about 3 km. Hence the 
Chapman function Is about 115 and the minimum detectable emission rate is about 
0.15 cm”*s *. This limit Is marked on all of the emission rate maps .and is 
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reflected in the iso-density profile maps for {O] as well. Another way of for- 
mulating the criterion for a minimum detectable emission rate makes use of the 
fact that the photometer's field extends 5 km in the vertical near 90 km. 

Hence* unless the emission rate averaged over 5 km is above a certain minimum 
value it cannot contribute significantly to the observed slant emission rate. 
That is to say that a condition that must be assigned is 

2/?Rd ~Y~ » Airl 

m Li 

where d is 5 km. R is the radius of the earth plus 90 km and 4irl is the minx- 

1ft 

mum detectable slant emission rate near 90 km. If the minimum detectable 
signal is* say 25R* because of uncertainty in the backgrouiid subtr'''tlon 

' E~ » 0.5 cm“’s”* 

D 

A glance at the emission rate viaps is enough to demonstrate that there is 
a large, amount of structure in the emission, both in altitude and latitude. 

The altitude of maximum emission varies with latitude but stays usually be- 
tween 95 km and 100 km. Most striking is the variation in latitude of the 
vertical gradient in emission rate. Regions in which the emission rate in- 
creases rapidly with altitude near 90 km alternace with regions of much 
slower variation. Vhe latitude difference between regions of rapidly vary- 
ing emission rates and slowly varying rates is 6* to 7*. In most cases where 
the emission changes slcwly, the emission rate deduced at 80 km is still sig- 
nificantly large by any criterion. In Figures lla-18a we show the resulti' of 
an error analysis of the calculated emission rates in addition to indicating 
on the graphs minimum observatlonally meaningful levels of emission rates: 

0.35 and 0.5 cm“*s”*. It can be seen that In the regions of steep gradients 



the emission rates increase rapidly near 90 km from values less than 0,5 cm" 
s”* to roaxlraums near f 7 km of the order of 20 to 100 On the other 

hand, at the latitudes where the emission rate varies slowly it often is in 
the range of 5 to 10 at 80 kn. Furthermore, the calculated values 

in these latter cases are usually physically significant according to the 
results of the error analysis. Pass 2269 (Figure 17) is particularly 
interesting in this regard. (This is the case in which 300 Iterations were 
performed.) In Figure 17a the emission rate profiles at 16.6* and 28* vary 
slowly from about 6 to 7 cm”*s“' at 80 km to about 40 cm-’s”* between 95 
and 100 km. On the other hand, at 18.6* and 23.6* the emission rates cal- 
culated at 88 km and below are mostly in the range 10”^ to 10'"*cm'"*s“* , far 
too small to be meaningful. By 90 km, they jump to about 5 cm"*s”*. 
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DISCUSSION 


The variation in the oxygen density with altitude and latitude will resemble 
the emission rate variation, the exact values Being determined by the mechanism 
selected for the green line excitation. The data for pass 2269 between 13* and 
28* were obtained at a local time of 24\cP, In Figures 17b, c, and d we show 
the densities at the 4 selected latitudes obtained following the analysis of 
Donahue et al, 119733. For comparison we also show the results of a calculation 
by S. C. Liu Iprivate conmmnication, 1977J using a time'<^ependent one~dimensional 
coupled flow model for atomic oxygen at local midnight and at 16Q0 hours. The 
time dependent equations that determine the behavior of odd oxygen between 8 Q 
and 120 km are the following: 


3I0J a^Oj,) 

- -[0]{ k^^[0H] + k^oCHO^l + k^glNO^l 


+ 2 k3g[0][M] + 2 k 35 [ 03 l) - [O^H k^^^lH] + k^-glOH] 


( 12 ) 


+ k23[H023> + 2 Jj^[02l 


3 [ 03 ] 3(J)(C3) 


3t 


3z 


= k3^[0)[02][M] - [03){k3^[H] + k33[0] 


+ kj^g[0H]} - 


(13) 


= -IO]{k3^[02][M] + 2 k3g[0][M]> + J2[03] + 2 J3[02] (14) 
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where the rate constants are defined in the paper by Liu et al . [1976] and 
the values used for these rate constants are those recommended by Hampson et 
al , [1977] and Mar?>itan ,^1977] . In this model, the solar zenith angle at noon 
was 25". The eddy diffusion coefficient above 80 km was held at 5.6 x lO^cm* 
s”*. The large diurnal variation In atomic oxygen below 90 km is caused by 
the fact that the conversion of 0 to 0^ by three body recombination 

0 + 02 + M->-02 + M 
is not compensated for by photolysis 

hV + 0^ 0 + O 2 

(The effect of 0 ^ photolysis at these altitudes is relatively Insignificant) . 
However, because of the non-negligible value for the eddy diffusion coefficient 
there is a strong downward flow of atomic oxygen below 90 km that causes the 
oxygen densities to be larger at any given altitude than the value that would 
prevail if photochemistry alone had determined the densities. Chemical loss 
rates for 0, 0^ and 0^^ as well as production rates are shown in Tables la 
and lb. These are gross rates to be used in a continuity equation 

- div (j> = p - 1 

where -is the oxygen flux, p the rate of 0 production and 1 the rate of 0 
loss. 
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Table la. Gross Production per cm’ per day 


Altitude 


Ox 


“3 

0 


110 km 

1.71 

X 10*® 

1.05 

X l(f 

1.71 X 

10*® 

110 

2.10 

X 10*® 

1.64 

X 10* 

2.01 X 

10*® 

90 

5.54 

X 10*® 

3.27 

X 10*® 

4.19 X 

10*® 

80 

1.00 

X 10** 

3.76 

X 10*® 

6.89 X 

10*0 

70 

2.79 

X 10?* 

1.67 

X 10* » 

1.15 X 

10*0 



Table lb. Gross 

Loss per cm® per day 


Altitude 


0 

X 


“3 

0 

110 km 

K74 

X 10^ 

1.04 

X l(f 

1.69 X 10* 

100 

1.01 

xlC*® 

1.64 

X 10* 

9.31 X 10* 

90 

6.96 

X 10** 

3.27 

X 10*" 

5.61 X 10*® 

80 

1.06 

X 10** 

3.90 

X 10*® 

7.58 X 10*® 

70 

2.78 

X 10** 

1.67 

X 10** 

1.15 X 10** 


It is very interesting that a characteristic of the derived densities in 
the latitude regimes of steep gradients is that they increase more rapidly with 
altitude between 85 and 90 km than the values computed from the coupled flow- 
chemical theory and much more slowly in the other regimes. This suggests that 
the effectiveness of downward transport varies strongly with latitude. At 18.6* 
and 23.6” on pass 2269 there appears to be weak downward transport of oxygen be- 
low 95 km. The midnight oxygen densities correspond to photochemical values. 

On the other hand at 16.6® and 28* during that orbit the downward transport 
is strong enough to cause midnight oxygen densities at 80 km orders of magnitude 




larger than those ext>ected from the ttoe dependent coupled flow model and even 
about 4 times as large as the daytime values predicted by theory. The simplest 
explanation we can find for these results is that downward transport in the 
neighborhood of the mesopause varies on a 6* latitude scale from very weak 
to very strong. 

One might have expected to find a reflection of this pattern in the magni- 
tude of the maximum 0 density. Where transport is strong, the maxinnim density 
should be relatively small. This expectation is disappointed. However, there 
is no reason to believe that a given flow condition is constant with altitude 
or that it is confined to a vertical surface. It may well vary with altitude, 
latitude and time. . It is probably significant that the structure in the slant 
emission rate as a function of altitude [Donahue et al ., 1973] shows a varia- 
tion on a scale of the order of 15*, as does the bottom side structure reported 
here. . 

• t 

In Figures 18 and 19 we show the results obtained On one Interesting pass, 
2505, on day 329 relatively early in the evening (2120 hours) and at high 
latitude. The very slow altitude variation between 80 and 104 km at 53.8* 
latitude is not characteristic of results obtained in any other case studied. 

A comprehensive summary of the maximum atomic oxygen densities (for case 
for day 304 and 311, as well as others, was presented by Donahue et al . 
[1974]. These data are given in the form of Iso-denslty contour maps. The 
nature of these maps indicates very large variability both temporally and 
spatially. Orbits which appear to overlap latitude regimes in this work, 
encompass regions which can be quite different as is seen in the patchy 
structure of the maximum atomic oxygen densities presented in the earlier 
paper. In the case of day 304 orbit 2142 generally has larger densities than 
orbit 2135 but is separated by about 180* of longitude and about 11 hours GMT, 
For the situation on day 311 of orbit 2238 and 2239 little change had taken 
place (see Donahue et al .. [19741) and the maps shown in Figure 4 and 5 are 



very similar in structure and magnitude. A dramatic example is the closed 
contour of about 4.8 x 10 atoms/r.m occurring at 20* latitude and 110 km 
in both Figures Aa and 5a. The consecutive orbits 2243 and 2244, however, 

Donahue et al . (1974] showed a very definite transition in densities. This 
can be seen in Figures 6 and 7 as decreased densities at all altitudes above 
97 km from orbit 2243 to 2244. The temporal variation cannot be resolved 
from the spatial in this analysis, but used in conjunction with the previous 
works (Donahue et al . [1973, 1974]) provide a good survey of the four 
dimensional variations of atomic oxygen. 
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Figure Captions 


Figure 1 . Geometry of the airglow observations. The outermost circle represents 
the satellite orbit, the next the green airglow emission maximum and the third 
the earth's surface. The line of sight of the photometer for selected scanning 
mirror positions is indicated at four successive positions of the satellite at 
which the green line is being observed. 

Figure 2a . Iso-emission rate contours in photons cmT^s'"* obtained during orbit 
2135 of OGO-6 on day 304, 1969. The heavy line at the top of the figure defines 
the transition to a region where the emission rate is lower than 0.15 cm~’’s~‘^ 
and the heavy line at the bottom the altitude at which the average emission 
rate becomes less than 0.5 cm~^s'~*. 

% 

Figure 2b . (Cj^) Contours of constant oxygen density on orbit 2135 in units of 
lO^^^cm"^ calculated according to the version of Che Qmpman mechanism used by 
Donahue et al . [1973]. Heavy lines separate regions tdiere the densities are 
calculated from meaningful and non-significant emission rates. 

Figure 2c . (C 2 ) Contours of constant oxygen density on orbit 2135 in units of 
10*®cm"* calculated according to the version of the Chapman mechanism developed 
by Wasser et al . [1976]. Heavy lines separate regions where the densities are 
calculated from meaningful and non-significant emission rates. 

Figure 2d . (SB) Contours of constant oxygen density on orbit 2135 in units of 
10 ’ “cm"’ calculated according to the version of the Barth mechanism developed 
by Slanger and Black [1977] . Heavy lines separate regions where the densities 
are calculated from meaningful and non-significant emission rates. 


Figure 3a. See Figure 2a» Iso-emlsslon rate contours In units of photons erne’s"*, 
orbit 2142. day 304, 1969. 

Figure 3b . (cp See Figure 2b, Iso-denslty contours In units of 10*®cm“*, orbit 
2142, day 304, 1969, Chapman mechanism [ Donahue et al ., 1973]. 

Figure 3c . (cp See Figure 2c, Iso-denslty contours in units of 10* “cm”*, orbit 
2142, day 304, 1969, Chapman mechanism [ Wasser et al ., 1976]. 

Figure 3d . (SB) See Figure 2d, iso-denslty contours in units of 10*®cm”*, orbit 
2142, day 304, 1969, Barth mechanism [ Slanger and Black , 1977). 

Figure 4a . See Figure 2a, iso-emission rate contours in units of photons cm~^ 

8“*, orbit 2238, day 511, 1969. 

Figure 4b . (C^^) See Figure 2b, iso-denslty contours in units of 10*®cm"*®, orbit 
2238, day 311, 1969, Chapman mechanism [ Donahue et al ., 1973], 

Figure 4c . (C2) See Figure 2c, iso-denslty contours in units of 10*®cm~’, orbit 
2238, day 311, 1969, Chapman mechanism, [ Wasser et al ., 1976). 

Figure 4d . (SB) See Figure 2d, iso-denslty contours In units of I0*®cm“*, orbit 

2238, day 311, 1969, Barth mechanism [ Slanger and Black , 1977). 

Figure 5a . See Figure 2a, iso-emlsslon rate contours in units of photons cm"^ 
s”*, orbit 2239, day 311, 1969. 

Figure 5b . (Cj^) See Figure 2b, Iso-emisslon contours In units of 10*®cm‘’*, orbit • 

2239, day 311, 1969, Chapman mechanism [Donahue et al., 1973). 


tiiiiiiiiffiiiiiMiiiisiiiN 


Figure 5c » (C 2 ) See •Figure 2c, iso^eneity contours in units of 10*®car*, orbit 
2239, day 311, 1969, Chapman mechanism [ Wasser et al. , 1976]. 

Figure 5d. (SB) See Figure 2d, iso-density contours in units of 10*®cbT®, orbit 
2239, day 311, 1969, Barth mechanism [ Slanger and Black, 1977]. 

Figure 6a. See Fi^re 2a, iso-emission rate contours in units of photons cmT^s"*, 
orbit 2243, day 311, 1969. 

Figure 6b. (C^) See Figure 2b, iso-density contours in units of 10^®cm*~®, orbit 
2243, day. 311, 1969, Chapman mechanism [ Itonahue et al. , 1973]. 

Figure 6c . (C 2 ) See Figure 2c, iso-denaity contours in units of 10*®cmT'*, orbit 
2243, day 311, 1969, Chapman mechanism I Wasaer et al ., 1976]. 

Figure 6d. (SB) See Figure 2d, iso-density contours in units of 10*®cm*'®, orbit 

2243, day 311, 1969, Barth mech?.nism [ Slahger and Black , 1977]. 

Figure 7a . See Figure 2a, iso-et:'.sslon rate contours in units of photons cm“*s“‘, 
orbit 2244, day 311, 1969. 

Figure 7b. (cp See Figure 2b, iso-denslLy contours in units of 10*®cm”*, orbit 

2244, day 311, 1969, Chapman mechanism [ Donahue et al . , 1973]. 

Figure 7c . (C 2 ) See Figure 2c, iso-density contours in units of 10'®cm“®, orbit 
2244, day 311, 1969, Chapman mechanism [ Wasser et al ., 1976]. 

Figure 7d . (SB) See Figure 2d, iso-density contours in units of 10'®cm“®, orbit 
2244, day 311, 1969, Barth mechanism [Slanger and Black, 1977]. 


Figure 8a . See Figure 2a, iso-eniaslon contours in units of photons caT^8~^, 
orbit 2269, day 311, 1969. 

Figure 8b . (C^) See Figure 2b, iso-density contours in units of 10* orbit 
2269, day 311, 1969, Chapman mechanism [ Donahue et al ., 1973]. 

Figure 8c . (Cj) See Figure 2c, iso-density contours in units of 10*®cm“*, orbit 
2269, day 311, 1969, Chapman mechanism [ Wasser et al ., 1976]. 

Figure 8d . (SB) See Figure 2d, iso-density contours in units of 10*®cm”*, orbit 
2269, day 311, 1969, Barth mechanism [ Slanger and Black , 1977]. 

Figure 9a . See Figure 2a, iso-emission contours in units of photons cm“*s“*, 
orbit 2505, day 329. 1969. 

Figure 9b . (Cj^) See Figure 2b, iso-density contours in units of 10*®cm“®, orbit 
2505, da> 329, 1969, Chapman mechanism [ Donahue et al ., 1973]. 

Figure 9c . (C 2 ) See Figure 2c, iso-density contours in units of 10*®cm“*, orbit 
2505, day 329, 1969, Chapman mechanism [ Wasser et al ., 1976]. 

Figure 9d . (SB) See Figure 2d, iso-density contours in units of 10*®cm”*, orbit 
2505, day 329, 1969, Barth mechanism [ Slanger and Black , 1977]. 

Figure 10 . &nlsslon rates in photons cm”*s“* for the same data as Figure 8a 
except that 300 Iterations were performed in arriving at this emission matrix 
Instead of 25 as in Figure 8a. 


Figure 11a . 557.7 nm volume emission rate in photons cb‘*’s‘~* as a function of 
altitude at selected latitudes of orbit 2135, day 304, 1969. Minimum detectable 
enission rates of 0.15 cm'^s**! and 0.5 cm~>s*’i are indicated by vertical lines. 
Error bars are assigned to calculated values at -26.6* latitude (triangles). 

Figure 11b . (C^) Oxygen densities calculated from the emission rates of Figure 
11a following the Chapman mechanism according to Ik>nahue et al. [1973]. Minimum 
meaningful densities corresponding to the minimum detectable emission rates of 
Figure 11a are indicated by vertical lines. 

Figure 11c . (C 2 ) Oxygen densities calculated from the emission rates of Figure 
11a folloid.ng the Chapman mechanism according to Wasser et al. [1976]. Minimum 
meaningful densities corresponding to the minimum detectable emission rates of 
Figure 11a are indicated by vertical lines. 

Figure lid . (SB) Oxygen densities calculated from the oDlssion rates of Figure 
11a following the Barth mechanism according to Slanger and Black [1977]. Mini- 
mum meaningful densities corresponding to the minimum detectable onission rates 
of Figure 11a are indicated by vertical lines. 

Figure 12a . See Figure 11a, orbit 2142, day 304, 1969. 

Figure 12b . (C^^) See Figure 11b, orbit 2142, day 304, 1969, Chapman mechanism, 
fbonahue et al . , 1973]. 

Figure 12c. (C 2 ) See Figure 11c, orbit 2142, day 304, 1969, Chapman mechanirai, 
[Wasser et al., 1976). 


Figure 12d . (SB) See Figure lid, orbit 2142, day 305, 1969, Barth mechanism, 
iSlanger and Black , 1977). 

Figure 13a . See Figure 11a, orbit 2238, day 311, 1969. 

Figure 13b . (C^) See Fi^re 11b, orbit 2238, day 311, 1969, Chapman mechanian, 

( Donahue et al. , 1973]. 

Figure 13c . (Cj) See Figure 11c, orbit 2238, day 311, 1969, Chapman mechanism, 
( Wasser et al ., 1976). 

Figure 13d . (SB) See Figure lid, orbit 2238, day 311, 1969, Barth mechanism, 
( Slanger and Black , 1977). 

Figure 14a . See Figure 11a, orbit 2239, day 311, 1969. 

Figure 14b . (C^ See Figure 11b, orbit 2239, day 311, 1969, Chapman mechanism, 

( Donahue et al . , 1973] . 

Figure 14c . (Cj) See Figure 11c, orbit 2239, day 311, 1969, Chapman mechanism, 
( Wasser et al ., 1976]. 

Figure 14d, (SB) See Figure lid, orbit 2239, day 311, 1969, Barth mechanism, 
( Slanger and Black , 1977]. 

Figure 15a. See Figure 11a, orbit 2243, day 311, 1969. 


Figure 15b . (Cj^) See Figure lib, orbit 2243, day 311, 1969, Chapman mechanism, 
(Donahue et al., 1973]. 


Figure 15c. (C 2 ) See 'Figure 11c, orbit 22A3, day 311, 1969, Chapman nechanism, 
I Wasser et al., 1976]. 

Figure 15d . (SB) See Figure lid, orbit 2243, day 311, 1969, Barth mechanisa, 
[ Slanger at^ Black , 1977). 

Figure 16a . See Figure 11a, orbit 2244, day 311, 1969. 

Figure 16b . (C^ See Figure lib, orbit 2244, day 311, 1969, Chapman mechanism, 

[ Donahue et al. , 1973]. 

Figure 16c. (C 2 ) See Figure 11c, orbit 2244, day 311, 1969, Chainaan mechanism, 
[ Wasser et al , 1976) . 

% 

Figure I6d . (SB) See Figure lid, orbit 2244, day 311, 1969, Barth mechanism, 
[ Slanger and Black, 1977). 

Figure 17a . See Figure 11a, orbit 2269, day 311, 1969. 

Figure 17b. (Cj^) See Figure lib, orbit 2269, day 311, 1969, Chapman mechanism, 

[ Donahue et al .. 1973). Model densities are shown for 1600 hours and 2400 hours. 

Figure 17c . (C 2 ) See Figure 11c, orbit 2269, day 311, 1969, Chapman mechanism, 

[ Wasser et al.. 1976). Model densities are shown for 1600 hours and 2400 hours. 

Figure 17d . (SB) See Figure lid, orbit 2269, day 311, 1969, Barth mechanism, 
[ Slanger and Black , 1977] . Model densities are shown for 1600 hours and 2400 


hours . 


Figure 18a . See Figure 11a, orbit 2505, day 329, 1969. 


Figure 18b . (C^) See Figure 11b, orbit 2505, day 329, 1969, ChaiMaan mechanism, 
t Ponahue et al .» 1973]. 

Figure 18c . (C^) See Figure 11c, orbit 2505, day 329, 1969, Chapman mechanism, 
[Wasser et al ., 1976]. 

Figure 18d . (SB) See Figure lid, orbit 2505, day 329, 1969, Barth mechanism, 
[Slanger and Black , 1977]. 

Figure 19 . Oxygen densities of Figure 18b compared to calculated values at 12^00 


16^00 and 20^)0. 
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